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Universal properties of the statistics of stochastic efficiency for mesoscopic time-reversal symme¬ 
try broken energy transducers are revealed in the Gaussian approximation. We also discuss how 
the second law of thermodynamics restricts the statistics of stochastic efficiency. The tight-coupling 
(reversible) limit becomes unfavorable, characterized by an infinitely broad distribution of efficiency 
at all times, when time-reversal symmetry breaking leads to an asymmetric Onsager response ma¬ 
trix. The underlying physics is demonstrated through the integer quantum Hall effect and further 
elaborated in a triple-quantunr-dot three-terminal thermoelectric engine. 
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Introduction .— Nonequilibrium phenomena of macro¬ 
scopic systems are described by their responses to exter¬ 
nal perturbations and follow the laws of thermodynam¬ 
ics [l|. Statistical fluctuations of a measurable quantity 
are negligible in macroscopic systems but in very small 
(“mesoscopic”) systems they play an essential role [i|. 
Importantly, in such a mesoscopic regime a single mea¬ 
surement of a physical quantity (e.g., electrical current) 
does not follow the laws of thermodynamics (although its 
average over many different measurements does). Well 
known examples are the Jarzynski equality and 

the fluctuation theorem, which states that a stochastic 
negative entropy production may show up, though expo¬ 
nentially unlikely compared to its corresponding positive 
entropy production process US- An important conse¬ 
quence of the fluctuation theorem is that the “stochastic 
energy efficiency” of a very small engine can be larger 
than the Carnot efficiency, although the average effi¬ 
ciency over many measurements is smaller or equal to the 
Carnot efficiency, in accord with the second law of ther¬ 
modynamics [il.l2.H0l. 12, 14 -171. The study of stochastic 
efficiency of very small energy transducers is pertinent to 
understanding of mesoscopic thermoelectric energy con¬ 
version in normal electron Hill and Cooper pair is¬ 
lands Hi biological photosynthesis in an individual re¬ 


21] , and photo-mechanical energy conversion 


action unit 
in quantum optomechanical systems 


Analyzing stochastic efficiency, it was recently shown 
that the Carnot efficiency is the least likely stochastic 
efficiency 0. later found to be solely the consequence 
of the fluctuation theorem Hi f° r time-reversal symmet¬ 
ric (TRS) energy transducers [1]. Breaking time-reversal 
symmetry can shift the least likely efficiency away from 
the Carnot efficiency H, H- However, little is known 
about efficiency statistics of time-reversal symmetry bro¬ 
ken (TRB) mesoscopic energy transducers, which is the 
gap we want to fill in this Letter. Another fundamental 
question concerns how the second law of thermodynamics 


restricts the statistics of efficiency for both TRS and TRB 
energy transducers. In this Letter we address these prob¬ 
lems for systems operating in the linear-response regime 
where fluctuations can be well described within the Gaus¬ 
sian approximation. One of our key results is that in 
the reversible (or “tight-coupling”) limit the distribution 
of stochastic efficiencies becomes infinitely broad at all 
times , thus, the average efficiency loses its meaning even 
in the infinite long time limit. This anomaly occurs only 
for reversible TRB energy transducers of which the On¬ 
sager response matrix is asymmetric. We discuss this 
anomaly using the example of the integer quantum Hall 
effect. We further demonstrate our result with a triple¬ 
quantum-dot (QD) thermoelectric engine where a full¬ 
counting statistics method confirms our Gaussian theory. 


Efficiency statistics in the Gaussian approximation .— 
We consider a generic situation in which there are two en¬ 
ergy output channels (“1” and ”2”). Each of the channels 
has a thermodynamic “current” and a conjugated driv¬ 
ing force (i.e., affinity). The time-integrated currents are 
denoted by Ji (i = 1,2) while the time-intensive current 
is defined as 7, = Ji/t with t the total time of opera¬ 
tion. The affinities are associated with the properties 
of the reservoirs (e.g., temperatures and electrochemi¬ 
cal potentials) and hence fluctuate negligibly. In con¬ 
trast, the currents may fluctuate considerably. A small 
TRB machine can be characterized in the linear-response 

regime by Ii = MijAj ( i,j = 1,2), or I = MA with 


I = (7i, 1 2 ) and A = (Ai,A 2 ). As well, in this regime 
the statistics of the currents at long time t can be de¬ 
scribed within the Gaussian approximation by the distri¬ 
bution P t (7) = ‘ Vdet(( f; 1)l ^ exp(—1(57^ • M~ l • 61) 
241 ]. Here det((A7 -1 ) sym ) is the determinant of the 


symmetric part of the inverse of the Onsager response 
matrix M and the superscript “T” denotes transpose. 
While averaged quantities are represented with a bar over 
the symbols throughout this paper, 61 = 7 — 7 repre- 
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sents fluctuations of the currents. From the probabil¬ 
ity distribution of stochastic currents we calculate the 
distribution of efficiency Pt{q) [!]• We then obtain the 
large deviation function (LDF) of the stochastic efficiency 
Q{q) = — limt-Kxj t~ l ln[P t (? 7 )]. The scaled LDF is de¬ 
rived in the Supplementary Material. It is given by 


m $ ft 

df.ot 


J{qc) (q + a 2 + aqr + aqq)' 


(1 + a 2 + aqr + aq) {rj 2 + a 2 + aqr] + aqrrj) 


, (1) 


where Stot = 'Yln^iAi I s the average total entropy pro¬ 
duction rate and 


J(vc) = 


4- q 2 { 1 + r) 2 
16(1 — q 2 r) 


( 2 ) 


is the scaled LDF at Carnot efficiency. Here, 


Mo 


q = 


\J M22M11 


M v 

M- 


a 


: 


A\\JM\\ 
A2 \J M22 


( 3 ) 


are dimensionless parameters that characterize the re¬ 
sponses of the system and the applied affinities. We term 
q the degree of coupling S3 , r the TRB parameter 0 , 
and a the affinity parameter Q. In addition, the effi¬ 
ciency is defined as rj = — IiAi/{l 2 A 2 ) MSI El. For 
thermal engines, q = fj/rjc , with the standard definition 
of energy efficiency rj = W/Q and rjc the original Carnot 
efficiency. In our scheme, efficiency is scaled so that the 
Carnot (reversible) efficiency corresponds to rjc = L 
The second law of thermodynamics requires that 
Stot > 0 and hence Mu, M 22 > 0 and M 11 M 22 > (M 21 - 
Mi 2 ) 2 /4, i.e., 


— 2 < q{l + r) < 2. (4) 

Equality is attained only in the “tight-coupling” limit EH 
where the average efficiency reaches its upper bound. It 
has been proposed El that breaking time-reversal sym¬ 
metry can open the possibility of achieving Carnot effi¬ 
ciency at finite output power and improving the efficiency 
at maximum output power to overcome the Curzon- 
Ahlborn limit S3- hr contrast to the average efficiency 
discussed in previous works 1, @,! Si, Si S3, the 
LDF J{rj) allows us to examine the statistics of efficiency 
fluctuations. Particularly, we will show that in TRB sys¬ 
tems the tight-coupling limit becomes unfavourable as 
the efficiency distribution becomes infinitely broad at all 
times. 

The LDF in TRB systems, Eq. |T|), is a key expres¬ 
sion in our work. As we discuss below, its shape can 
be characterized by three quantities: the average value 
of efficiency rj, the least probable efficiency q*, and the 
width of the distribution around the average, a v . In what 
follows we investigate these properties, particularly, un¬ 
der two (separate) experimentally relevant conditions of 



FIG. 1. Efficiency statistics for TRS and TRB systems, (a) 
An example of the LDF J(q) for TRS (blue, r = 1) and TRB 
(green, r = 1.6) cases with a = —0.3 and q = 0.7. (b) The 
least probable efficiency q* for these two systems at different 
a. (c) LDF J(i]) as a function of a and q for a TRB system 
with r = 1.6 and q = 0.7. Note that the least probable 
efficiency (the maximum of the LDF) q* shifts away from the 
Carnot efficiency [qc = 1, labeled as the dashed line in (c)] 
for TRB systems as demonstrated in the panels (a), (b), and 
(c). In panel (c) the minimum of J(q ) for q 1 with large 
negative a is associated with the reversed machine of which 
the actual efficiency is l/q, instead of q. (d) Our analysis is 
exemplified on a triple-QD thermoelectric device, see Fig. [4] 
and text. 


maximum average efficiency and maximum average out¬ 
put power. 

We begin with some general properties of the LDF. 
First, the LDF at the Carnot efficiency J{qc ) is inde¬ 
pendent of affinities but it is solely determined by the 
response coefficients. It is also invariant under time- 
reversal operation, which turns r —> 1/r and q —> qr 
S3, 261. J{qc ) can be suppressed by breaking time- 


reversal symmetry, particularly in approaching the tight- 
coupling limit. Second, in the TRS limit, r = 1, Eq . (HD 
goes back to results obtained in Refs, [l], S, 0, 0. In 
more general situations we find that 0 < J{q ) < 1/4 
is guaranteed by the thermodynamic bound d4} [see Sup¬ 
plementary Material]. Moreover, J{q) has only one mini¬ 
mum and one maximum. While the minimum J{rf) = 0 is 
reached at the average efficiency rj = -a(a+qr)/(aq- 1-1), 
the maximum value J(q *) = 1/4 is realized at the least 
probable efficiency 


1 + 


q{i — 1)(1 + aq + aqr + a 2 ) 
q — qr — 2a + q 2 { 1 + r) 


( 5 ) 


In the TRS limit, the least likely efficiency is always iden- 
tical to the Carnot efficiency, q* = qc = 1 MlHEI. 
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For TRB systems, in contrast, we find here that 77 * de¬ 
pends on the parameters q , r, and a, see Fig. [2a), (b), 
and (c). We note that rj* diverges at a c = 2 -<^q* r • 
For |r| > 1, a c produces a positive average efficiency 
and output power, relevant for device operation. Be¬ 
sides, 77 * r ~ > °°> oo for all a and q. It has been shown 
that this limit is achievable for a triple-QD thermoelec¬ 
tric device [see Fig.QJd) for schematic of the devicelwhen 
M 21 —> 0 and M 12 ^ 0 take place simultaneously | 6 |, see 
also Supplementary Materials. 

The width of the distribution around the average effi¬ 
ciency, tr v , is another key characteristic of efficiency fluc¬ 
tuations. Expanding J{rf) around its minimum 77 , one 
writes J(rj) ~ 53^(77 — 17) 2 + 0((ry — 7 ?) 3 ), to provide here 


2 v / 2 |a|(l — q 2 r)( 1 + a 2 + aq + aqr) 
(1 + aq ) 2 ^/4 - q 2 ( 1 + r) 2 


( 6 ) 


We now proceed to describe the properties of 77 , 77 * and 
a v under conditions for optimized (average) operations. 

Efficiency fluctuations at maximum average 
efficiency. — We obtain the familiar form for the 
maximum average efficiency 


VZT + 1-1 

V = Vmax = r , - = r 


1- v/T - ? 


yj ZT +1 + 1 \ 1 + yj\ — q 2 r 


(7) 


when a = — qr/{\ + yjl — q 2 r ) Q, expressed in terms of 
the figure of merit for energy conversion ZT = q 2 r/(l — 
q 2 r) p, [2f|. The thermodynamic upper bound of the 
average efficiency, reached at the tight-coupling limit, 

| g( l + r)|^2,isiM 


Abound = min{r 2 , 1 }. 


( 8 ) 


One of our key results here is that under the maximum 
average efficiency condition, the least probable efficiency 
© reduces to ( 77 c = 1) 


77 * = r. (9) 

In the TRS limit this recovers recent findings that Carnot 
efficiency is the least probable stochastic efficiency 00 . 
The maximum average efficiency, the upper bound of the 
average efficiency, and the least probable efficiency are 
plotted in Fig. ©a) as a function of r for q = 0.5, and we 
observe that 


V — Vbound — V maxi (10a) 

V ^ Vmax < Vboundi Vr < 0. (10b) 

The least probable efficiency 77 * coincides with rjb OU nd 
only in the TRS limit, r = 1, or for the TRB case with 
r = 0 . 

The width parameter o v is plotted in Fig. ©b) where 
the white region in the figure is forbidden by the second 




FIG. 2. Efficiency statistics under the maximum average effi¬ 
ciency condition, (a) Maximum average efficiency T 7 maa , (blue 
curve), thermodynamic upper bound of the average efficiency 
Vbound ( re( l curve), and the least probable efficiency rf (green 
curve) for q = 0.5. The blue dot represents the TRS point, 
r = 1. (b) Width of efficiency distribution Or,. The dashed 
line represents the TRS limit, r = 1. The white region is 
forbidden by the thermodynamic bound di- 


law of thermodynamics according to ©. It is small when 
q or r are small, corresponding to the “weak-coupling 
limit” when the average efficiency is small. Physically, 
this can be understood as that two nearly uncoupled 
currents are very unlikely to have collective fluctuations 
which result in a considerably large efficiency. 

Approaching the tight-coupling limit |g(l + r)| —> 2 
[which cannot be depicted in Fig. ©b)], we find that 
o v —> 00 for 0 < |r| < 1 and r = — 1 , while —> 0 
for other regimes. In the TRS limit, r = 1, our results 
agree with Ref. _l|. The singular behavior of o v in this 
limit can be understood by noticing that the denomina¬ 
tor of Eq. © vanishes in the tight-coupling limit whereas 
the numerator is proportional to the (average) total en¬ 
tropy production rate. It has been shown 0 [26| that 
in the tight-coupling limit, for |r| < 1 , the maximum 
average efficiency is attained at finite (average) total en¬ 
tropy production rate (thus the upper bound efficiency is 
less than 100%). In contrast, for |r| > 1, the maximum 
average efficiency is reached when the (average) total en¬ 
tropy production rate is zero (hence the upper bound 
efficiency is 100%). When S to t = 0 the unsealed LDF 
G(v) = Stot J ( 77 ) is always zero. Therefore, in the tight- 
coupling limit the distribution of stochastic efficiency is 
infinitely broad for |r| > 1 as well. The only exception 
is the TRS limit, r = 1, where the width of efficiency 
distribution becomes zero[l). 

Direct examination of Eq. m shows that in the tight- 
coupling limit J(jqc) = 0 (except for r = 1, the TRS 
limit). J{jf) thus vanishes whenever the average en¬ 
tropy production rate Stot is nonzero. However, when 
Stot = 0 , J(ji) is ill-defined, nevertheless the unsealed 
LDF Q(jj) = StotJ(v ) is always zero. In addition, for 
r = 0 and a = 0, J{rf) is constant for all 77 . Therefore, 
the distribution of efficiency is infinitely broad for any 
TRB energy transducer in the tight-coupling limit. 
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An example that may help to understand the anomaly 
in the tight-coupling limit is the integer quantum Hall 
effect. For example, at Filing factor u = 1, electrical 
transport in the quantum Hall insulator is described by 


e 

h 


0 1 
-1 0 


( 11 ) 


where j x / y and £ x / y are the electrical currents and fields 
along the x/y direction, respectively. Although the quan¬ 
tum Hall insulator does not conduct electron longitudi¬ 
nally, charge can be conducted via the Berry phase ef¬ 
fect, or more physically through the chiral edge states. 
In this way, the system converts electrical energy in the 
x direction W x = j x £x = jf£ x £y to electrical energy in 
the y direction W y = j y £ y = — ^£ x £ y (and vice versa). 
The macroscopic efficiency y = —W y /W x is always 100%, 
and the output power — W y is finite for nonzero (positive) 
f f 

c 'x i -'y • 

However, in the Gaussian description, the distribu¬ 
tions of the electrical currents j x and j y are singular be¬ 
cause the transport is completely dissipationless. This 
also leads to singular distributions of the output power 
and efficiency. Casting into our parameters, the integer 
quantum Hall systems have q = oo and r = — 1. Never¬ 
theless, in a finite system the longitudinal conductance is 
not vanishing (i.e., q is finite). A further examination of 
current noises in a Hall measurement needs more careful 
treatment 34[, which is beyond the scope of this work. 


We speculate that close to the tight-coupling limit, the 
Gaussian approximation is inadequate to describe fluctu¬ 
ations in the system. However, energy efficiency is much 
easier to measure in this system since the measurement 
of stochastic electric currents in mesoscopic systems is a 
rather mature technology [34). It is much easier to mea¬ 
sure stochastic efficiency distribution in this time-reversal 
symmetry broken system than in other known solid-state 
systems. 

Efficiency fluctuations at maximum average output 
power. — We now turn our attention to another highly 
pursued situation, the maximum average output power 
condition [17j arrived at a = —qr/2. Under this condi¬ 
tion the average efficiency is [B|, l26j 




rZT 


2 2 
q r z 


2(ZT + 2) 4 — 2q 2 r ’ 

with the thermodynamic upper bound [H. 261 


VboundiWmax) — ^ 


+ H 


( 12 ) 


(13) 


reached in the tight-coupling limit. Note that for r 2 > 1 
the efficiency at maximum power can be larger than 50% 
(the value of the Curzon-Ahlborn efficiency in the linear- 
response regime). The least likely efficiency is found to 
be 


r) = r 


4 — 3 q 2 r — q 2 r 2 
4 — 2 q 2 r — 2 q 2 r 2 


(14) 



FIG. 3. Efficiency statistics under the maximum average 
power condition, (a) Least likely efficiency y* as a func¬ 
tion of r and q. The dashed line represents the TRS limit 
r = 1. The white region outside is forbidden by the thermo¬ 
dynamic bound (j4j). In contrast, the white region inside the 
color graphics region is due to divergences of y* at special r-q 
lines given by r = (—1 ± \J\ + 8/q 2 )/2. (b) Least probable 
efficiency rf , average efficiency rj(W m ax) and its upper bound 
fjbound(Wmax) at q = 0.5 as a function of r. The dots identify 
values in the TRS limit, r = 1 . 


Calculations in Fig. [3] indicate that there are lines of sin¬ 
gularities for 77 *, varying with q and r. These singularity 
lines appear at r = (—1 ± 1 + 8 /g 2 )/ 2 , and they reach 

to the thermodynamic bound in the r-q plane at r = 1 
and q = ±1. The width of efficiency distribution un¬ 
der maximum average output power condition also shows 
the same singular behavior as that under the the maxi¬ 
mum average efficiency condition. In the Supplementary 
Material we explore J(y) under different situations and 
manifest its rich features. 

TRB thermoelectric transport in three-terminal 
systems. — We exemplify our analysis within a meso¬ 
scopic triple-QD thermoelectric device, see Fig. □Id). 
The affinities for a two-terminal thermoelectric device 
are the electrochemical potential A\ = iJ ' L eT ^ j R and the 
temperature difference A 2 = — ^r-, where e is the 

electronic charge, Tpp and pl.r denote the tempera¬ 
tures and chemical potentials in the left ( L ) and right ( R ) 
electronic reservoirs. In order to receive M 12 ^ M 21 , we 
introduce a third (probe P ) terminal and employ the con¬ 
straints that the average thermal and electrical currents 
flowing out of the probe terminal are zero @,[11. These 
conditions set the temperature Tp and chemical potential 
Pp in the probe. Each QD is coupled through elastic tun¬ 
neling to the nearby reservoir thus we employ the indices 
1/2/3 to identify the leads L/R/P , respectively. Hopping 
between QDs are affected by the magnetic flux $ pierc¬ 
ing through at the center with <f> = 27r$/$o ($0 is flux 
quantum). The system is described by the Hamiltonian 
P = Pqd T Plead T P-tun where Pqd = ^2/i— 1,2,3 ~f 

(te^/ 3 4+irfi + H.C.), Hlead = £^1,2,3 £* £kc\ k c ik , and 

Htun = I Zi,k vAc* + H.c.. 

This noninteracting model has been analyzed thor¬ 
oughly in Ref. [fj| using the Landauer-Biittiker approach, 
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E 2 /(k B T) 


E 2 /(k B T) 


FIG. 4. Triple-QD thermoelectric systems at the maximum 
average output power condition: (a) Least probable efficiency 
rf and (b) width of efficiency distribution a v . The QD con¬ 
nected to the probe terminal is set at E 3 = 2 . <f> = n/ 2 , 

r = 0.5, and t = —0.2. The white region in (a) depicts very 
large or very small (negative) r/* values which are not properly 
incorporated into the figure. 


to study transport properties and the average efficiency. 
We recall here that the transmission function is given by 
S(u>) = — i -MrG r (uj) where G r (w) = / 2)1 — H q d] 

is the retarded Green’s function of the quantum dots and 
the damping rate T = 2 n^ 2 k |V)fc| 2 5{uj — etk) is assumed 
to be a constant (independent of energy) for all three 
leads. Using <§> = tt/2, T = 0.5 and t = —0.2 (the equi¬ 
librium chemical potential is set at zero, the energy unit 
is fcflT), we calculate transport coefficients and substi¬ 
tute them into Eq. m to obtain the LDF of stochastic 
efficiency. We then calculate the least probable efficiency 
77 * and the width of distribution a v under the maximum 
average output power condition, see Fig. [4] We find that 
around Ei+E 2 = 0, \r/*\ becomes very large. The under¬ 
lying physics is that when E\ + E 2 is close to zero, the 
transport coefficient I-M 21 I can become very small while 
\M\ 2 \ is still finite, to yield a very large |r|[6]. Accord¬ 
ingly, the least probable efficiency I 77 * | becomes very large 
as we showed in Fig. [3jb) . In addition, at this special 
region the width cr v becomes very small, in accordance 
with Eq. ©, for a = —qr/2 (|gr| -C 1, q —> 0). These 
calculations, based on the Gaussian approximation, agree 
with a careful full-counting statistics analysis with vertex 
corrections Q, carried out in the linear-response regime 
[see Supplementary Material]. 

Future perspectives Our analysis lays the founda¬ 
tion for studies of efficiency statistics beyond the linear- 
response approximation, as well as to experimental works 
via, e.g., mesoscopic quantum Hall systems^ or TRB 
thermoelectric engines. 
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SUPPLEMENTARY MATERIALS 
I. LDF IN TRB SYSTEMS AND ITS PROPERTIES 


IA. Derivation of the LDF for TRB systems 

The LDF of efficiency fluctuations was derived in Ref. [lj, specific to TRS systems. Here we extend this work 
and obtain the LDF for TRB systems. We begin by introducing the probability distribution function (PDF) of the 
stochastic currents p (i = 1,2) 


t\J det((M _1 ) S j /m ) t 

Pt(I) = ^-exp(- i <5T r • AT 1 • SI). 


(15) 


Replacing the stochastic currents with stochastic entropy production rates S) = liAi [i = 1,2) at given affinities Ai 
and A 2 , one finds the PDF of the entropy production [lj, y] 


Pt(S 1 ,S 2 ) = 


t\/ det((C' ^sym) 


4:irSt 


■ exp 


— J—5S T ■ C~ l ■ SS 

'iStot 


(16) 


where SS = S — S, S = (Si,S 2 ) T is the stochastic entropy production, S = (Si,S 2 ) is the averaged (macroscopic) 
entropy production rate, and det((C' _1 ) S j /rra ) is the determinant of the symmetric part of the inverse of matrix C. The 
macroscopic total entropy production rate is S to t = Si + S 2 . Here, C-ij = AiMijAj / S to t ( i,j = 1,2). The PDF of 
stochastic efficiency is 


Pt(V) = J dS 1 dS 2 S(ri+^)P t (S ll S 2 )= J dS 2 \S 2 \P t (-r)S 2 ,S 2 ). 

A direct calculation yields an expression 


D , _ t^det((C l ) sym ) 

Pt{-r]S 2l S 2 ) = - -= -exp 

47 iStnt 


4 St, 


with the coefficients 


[a(y)S% + 2b{rj)S 2 + c ] 


= S t . 


a{y) = (C u + y(C 12 + C 21 ) + 7 ? 2 C 22 )/det(C), 
b(v) = {Si(C 12 + C 21 ) - 2 S 2 Cu + y[2C 22 S 1 - (C 12 + C 2l )S 2 ]}/[2 det(C)]. 


(17) 


(18) 


(19a) 

(19b) 


The coefficients Cy sum up to unity J2ij = 1- det((7) = Ci\C 22 — C\ 2 C 2 \ is the determinant of C. Note that 
C 12 7^ C 21 because of time-reversal symmetry breaking. The full probability distribution of the stochastic efficiency 
is now found to be 


Ptiv) = 


det((C 1 ) sym ) exp[—f,S tot /4] 

7T a(y) 


1 + V TTtS to th(y ) exp [tStoth (?y)]erf(y tS tot h(rj)) 


( 20 ) 


with erf (a;) = J * e * dt being the error function and 


h(ri) = 


~Hv) 


2S to t\/ a(y) 


( 21 ) 
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The large deviation function of stochastic efficiency is obtained from 

i t , w 

t^tot 4 

By substituting the following parametrization, 

M 2 i 


g = 


M12 Ai^/ Mu 

r = --, a = 

A / A A 

22 


a 2 Vm, 


we find that 


a (y) = 


M 2 i 

(1 + aq + arq + a 2 )(a 2 + aq( 1 + r)rj + rj 2 ) 
a 2 (l — q 2 r ) ’ 

2 (1 + aq + arq + a 2 )(a 2 q(r — 1) + a(q 2 ( 1 + r)(r - 77) + 2 (rj - 1 )) + q(r - l)r]) 


b{q) = M 2 2 A 2 - 

Stot = M 2 2 Al(l + aq + arq + a 2 ) 

J(v) = 


2 a(l — q 2 r) 

[4 — q 2 ( 1 + r) 2 ] ( a 2 + rj + aqr + aqr )) 2 


16 (1 — q 2 r) (1 + a 2 + aq + aqr) (a 2 + aqr] + aqrr] + rj 2 ) 


( 22 ) 

(23) 

(24a) 

(24b) 

(24c) 

(24d) 


IB. Thermodynamic bounds on the LDF 


We prove here that the inequalities 0 < J{y) <1/4 are guaranteed by the thermodynamic bound 

k(l + r)|<2. (25) 

First, 4(1 — q 2 r) > 4—q 2 (l+r ) 2 > 0, which guarantees the positive semi-definiteness of the prefactors of the numerator 
and denominator in Eq. (I24dl> . Second, 1 +a 2 +aq+aqr = (a + q A+A ^ -f ^1 — iSlALLjj > 0 . This is consistent with 
the fact that this term originates from the average total entropy production rate. The last term in the denominator 
is also not less than zero, since a 2 + aqrj + aqqr + r ] 2 = (a + + rj 2 ^1 — > (y Therefore J(rj) is 

guaranteed to be greater than zero. Using exactly the same arguments, one can show that 


I _ = [« 2 g(r - 1) + «(g 2 (l + r)(r -y)+ 2 (rj - 1 )) + q(r - l)r;] 2 > Q 

4 16 (1 — q 2 r) (1 + a 2 + aq + aqr) (a 2 + aqr] + aqrr] + rj 2 ) ~ 

Therefore, 0 < J(rj) < 1/4. From the above we find that J{rj) = 1/4 is reached only at 

q{i — 1)(1 + aq + aqr + a 2 ) 


ri = 1 


q — qr — 2 a + q 2 ( 1 + r)a 


(26) 


(27) 


We also find that there is only one minimum of Jjrj) which is the macroscopic efficiency rj , and only one maximum 
of J{rj) which is precisely the least probable efficiency t]* given above. This is confirmed by solving the extremum 
equation <9r;</(? 7 ) = 0 where we find only two solutions: one is rj , the other is rf. This property determines the 
basic-generic shape of the LDF curve. 


IC. J(q) under time-reversal operation 

In the main text we defined the following parameters 

lb/12 


M 2 i \/ M22M11 

Under time-reversal operation, cf> —» — </>, the above parameters transform as follows 

r (- < />) -t K g(- 0 )-tgr, 


(28) 


(29) 























where we have used Onsager’s reciprocity relation M^iA) = -W 21 (—</>)■ Denoting the LDF of efficiency for the reversed 
magnetic field {<j> —»• —<t>) by we obtain that 


J{iD - Av) 


J(rjc) a 2 q 2 {r 2 

1 + a 2 + aqr + aq 


1)(1 - rj 2 ) + 2 (a 2 + rj)aq{r - 1)(1 - 77 ) 
a 2 + aqr/ + aqrr] + rf 2 


(30) 


We find that for TRB systems Carnot efficiency (rjc = 1 ) appears as a special point where the distributions become 
invariant under time-reversal operation, i.e., J(j]c) = J(vc)- For TRS systems (r = 1) this equality trivially holds 
for all values of efficiency. 


ID. J(jj) at various limits 

In this section we illustrate the rich behavior of J{jf) at various limits: (i) weak coupling limit r —> 0 or q —> 0; (ii) 
tight-coupling limit with maximum macroscopic efficiency for 0 < |r| < 1 , |r| > 1 , and r = ±1; (iii) tight coupling 
limit with maximum macroscopic output power for |r| A 0,1, and for r = ±1. Results are plotted in Fig. [5] 

Fig. Eta) shows that in the weak-coupling limit (r —> 0), the LDF experiences a sharp transition from the minimum 
77 —> 0 to the maximum 77 * —> 0. Therefore, J{rj) behaves like a derivative of the Dirac delta function. In contrast, 
when q —^ 0 but r is finite, Fig. [5] (b) shows that J(rf) develops an infinitely narrow dip near 77 = 0, i.e., it behaves 
like the Dirac delta function itself. Note that in the weak coupling regime with q —> 0 and/or r —> 0, the affinity 
parameter a — > 0. As a result, the behavior of J/ 77 ) is very similar either under the maximum macroscopic efficiency 
condition or the maximum output power condition as both the efficiency and the output power go to 0 . 

In the tight-coupling limit, the behavior is quite different under those two conditions. We first examine the maximum 
macroscopic efficiency condition. Fig.EIc) shows that for 0 < |r| < 1 , the width of efficiency distribution a v tends to 
infinity, while the maximum value of J{jf) at 77 * = r develops a very sharp peak. In contrast, for |r| > 1, as shown in 
Fig. [He), the width of efficiency distribution approaches zero, while the width at the least probable efficiency becomes 
infinite. The r = 1 situation, see Fig. Eld), demonstrates a sharp transition from the minimum value to the maximum 
point, resembling the behavior of the derivative of the Dirac delta function. For r = —1, the tight coupling limit, 
|g(l + r)| — > 2 is pushed to q — > 00 . Therefore, in this situation for any finite q the behavior of J(jf) shows a regular 
behavior. Nevertheless, as shown in Fig. [5|f)j the distribution of 77 is quite broad for r = —1. This case is relevant 
to recent studies on “chiral thermoelectrics” (e.g., Nernst engines) where, however, a much stronger bound on q was 

obtained Hi- 

Under the maximum macroscopic output power condition, the behavior of J(rj) for r —> 0, q — > 0 cases is identical 
to that observed in Fig.[5|a) and (b). The tight-coupling limit with 0 < |r| < 1 also shows features similar to Fig. (He). 
However, when r = 1 in the tight-coupling limit, see Fig.EJg), the width around the macroscopic efficiency (77 = 0.5) 
approaches zero while the width around the maximum value of J{rf), at 77 * = 1, becomes infinite. Fig. (Hh) focuses on 
the |r| > 1 case for which the width of efficiency distribution tends to infinity, while the maximum of J(rf) becomes a 
sharp peak. FigEfi) shows that at r = — 1 , J/ 77 ) becomes an extremely broad distribution. 


IE. Width of efficiency distribution under general conditions 

The width of efficiency distribution at arbitrary a is found to be 

_ 2v / 2|a|(l - g 2 r)(l + a 2 + aq + aqr) 

^ ~ (1 + aq) 2 y/l^q 2 (l+r) 2 ' 

In Fig. [G] we display a v as a function of the affinity parameter a and the TRB parameter r for q = 0.5. We find that 
the width of efficiency distribution grows with |a|; at a ~ 0 it is very small. Besides, it tends to very large values in 
approaching the tight-coupling limit. 


IF. Width of efficiency distribution under the maximum macroscopic power condition 

Fig. 2(b) in the main text illustrates the behavior of the width of efficiency distribution <j v under the condition of 
maximum macroscopic efficiency. Here we complement this result and display in Fig. [7|a) the behavior of a v under 
the condition of maximum macroscopic power. We see in this figure that the features of <7 V are qualitatively similar in 
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FIG. 5. Illustrations of the rich behavior of J(ji) at various limits, (a)-(b) Weak-coupling limit, (a) The r —> 0 limit with 
r = 0.001 and q = 0.5. (b) The q —> 0 limit with finite r. We employ q = 0.001 and r = 0.5. (c)-(f) Several examples for 
achieving the tight-coupling limit |g(l + r)| —>- 2 with maximum macroscopic efficiency, (c) r = 0.5, q = 1.3332 (d) r = 1, 
q = 0.9999, (e) r = 5, q = 0.3332, (f) r = —1, q = 10. (g)-(i) Examples for the tight-coupling limit with maximum macroscopic 
output power, (g) r = 1, q = 0.9999, (h) r = 5, q = 0.3332, (i) r = —1, q = 10. Note that the scales of the horizontal axes in 
these figures are different. 



r 


FIG. 6. Width of efficiency distribution a v as a function of the affinity parameter a and the TRB parameter r for q = 0.5. The 
boundaries, r = —5 and r = 3, satisfy the tight-coupling limit. White regions appear either because a v is too small (around 
a = 0), or because it is too large (around r = —5 or 3). 


both cases. In particular, at the small q 2 r limit, <j v for both cases become quantitatively similar, which is consistent 
with the understanding that the maximum macroscopic efficiency and maximum macroscopic power conditions are 
close at this limit [H]. Nonetheless, we find that a v diverges in the tight-coupling limit for all r except at r = 0,±1. 
In Fig. 0(b) we examine a particular case with q = 0.5. The tight-coupling limit is then reached for r = —5 or r = 3, 
where a v is shown to diverge. 
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FIG. 7. width of efficiency distribution at the maximum macroscopic output power condition, (a) Width of efficiency 
distribution a v as a function of r and q. (b) Particular example of a v with r at q = 0.5. 




FIG. 8. Efficiency statistics while considering the quantum bound on transport coefficients, (a) Width of efficiency distribution 
a ,, under the maximum macroscopic efficiency condition, (b) Least probable efficiency rf and (c) width of distribution a v under 
the maximum macroscopic output power condition. 


II. QUANTUM MECHANICAL BOUND ON EFFICIENCY FLUCTUATIONS 


We examine here the effect of the quantum mechanical bound, which was introduced in Ref. Q, on efficiency 
fluctuations. Within our parameters the bound comes up as 1 + q 2 ?— q 2 — q 2 r 2 > 0, affecting the width of efficiency 
distribution fluctuation a v . The bound was obtained by considering a three-terminal Landauer-Biittiker model, and 
it is a direct result of the unitarity of the scattering matrix. In Fig. [ 8 ] we plot (a) <j v under the maximum macroscopic 
efficiency condition, (b) the least probable efficiency rf, and (c) under the maximum macroscopic power condition. 
The quantum bound regularizes the divergent behavior of cr v under both conditions, as well as the divergency of rf 
for the maximum power condition. This is because the above bound prohibits us from meeting the tight-coupling 
limit, unless r = 1, see J(rj) in Fig. [5jd). However, note that the least probable efficiency can still diverge at the 
following affinity parameter, even within the quantum mechanical bound, 


Oi c 


g(l - r) 

2 - q 2 - q 2 r ' 


(32) 


At this value we receive positive average efficiencies and average output power when r > 1 or r < - 1 . Inspecting 
Fig. [5] further, we find that the width of efficiency distribution under both maximum macroscopic efficiency and 
maximum macroscopic power conditions can become considerably large. In addition, the least probable efficiency rf 
can deviate significantly from the Carnot efficiency. In fact, within the maximum macroscopic efficiency condition, 
77 * is allowed to diverge when r —> 00 Q- We note that the above bound goes to the thermodynamic bound for 
TV-terminal Landauer-Biittiker conductors, when N —> 00 [!]. However, the quantum bound may not necessarily hold 
when genuine many body interactions and inelastic processes (e.g., electron-electron scattering and electron-phonon 
scattering) are taken into account. 


















III. COMPARISON BETWEEN THE GAUSSIAN APPROXIMATION AND A FULL COUNTING 
STATISTICS ANALYSIS FOR A TRIPLE-QDS THERMOELECTRIC MODEL 
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In this section we present a comparison between the LDFs as obtained under the Gaussian approximation and a 
full counting statistics analysis using a three-terminal thermoelectric model. A counting statistics method including a 
probe terminal has been recently developed by Utsurni et a 1. (though the discussion there is focused on heat currents 
fluctuations) Q- In this approach the LDF of efficiency is calculated through the cumulant-generating function 
(CGF) defining counting fields for the particle number and energy A j in each reservoir j = L,R,P. The CGF of 
the three-terminal system is given by [7| 

P 3 1 = J ^ lndet[i — <j>)]. (33) 


Here 

/ = diag(/ L ,/ p ,/ P ), K = 1 - e l§ S f {uj,(/))e~ l§ S(uj,(j)), 9 = diag(wA p + £ L , wA p + £ Ri ujX p + £ P ), (34) 

with “diag” denoting a diagonal matrix, fj = [exp( ^~ Mj ) + 1 ] _1 (j = L,R,P) is the Fermi distribution for j —th 

reservoir, and S is the S-matrix of the triple-quantum-dot (QD) system. 

By integrating out the short-time dynamics it was shown in Ref. Q that the effective CGF for particle and energy 
transport between the L and R reservoirs (when the particle and energy currents flowing out of the probe terminal 
are zero) is given by 

^ 2 t(AL, ^L,T R , fl R ) = J r 3t(\ L ,\ R ,\*p 1 t;Lit;RiCpiTL 1 H'L,T Rl l± R iTp,li*p). (35) 


Here X * P , Tp , and fip are determined by the saddle-point equations 

dP: it dP 3t dP :it dP u 

dX P d£, P dTp dfj,p [ ’ 

which maximize the probability of processes with zero energy and particle currents flowing out of the probe terminal. 

In the linear-response regime, the effective two-terminal CGF can be approximated by a second-order expansion in 
A j and £j as well as the affinities {jij — /J.)/T and 1/T — 1 / Tj for j = L,R. Due to particle and energy conservation, 
the counting fields for the right reservoir can be regarded as redundant, hence we can set X R = l; R = 0. Furthermore 
we set T r = T and p R = //. The approximate two-terminal CGF is now given by 

tip) = ^a T ■ Ha, (37) 

where 

d r = (A i ,a,A 1 ,A 2 ), H 80 , = 

with A 1 = (n L - h r )/T r and A 2 = 1/T R - l/T L . 

The matrix TZ must be calculated from the second derivative tensor of the three-terminal CGF at equilibrium with 
vertex corrections as shown in Ref. [7|, 



d 2 p 2 , 


dapdap' 


-(■ 


£> 2 Pii a 2 P 


3 1 


U~ 


d 2 T- 


3 1 


Xl—^L—Xp—^P—Ai—A2—A3—A4=0 


\ l =£ l =Ai=a 2 =o \dapdap' dapdbry db^da^ 
where b = (A P , £ P , A 3 , A 4 ) with A 3 = (/r P — fi R )/T R , A 4 = 1/T R — 1 /Tp, and U = B~ l with 

d 2 T 3t 


77 db-ydb~ 


(39) 


(40) 


Xl—^L — Xp=^P—Ai=A2—A3 — A4=0 


We now change variables 


iXl -t r]((T L - T r )/T 2 , i£ L -A £(ml - /r p )/T, 


(41) 
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FIG. 9. Comparison between J(ri) within the Gaussian approximation (labeled by “G”), and from a full counting statistics 
analysis (labeled by “FCS”) for the triple-QDs thermoelectric model, (a) J{rj) for cj> = n/2. (b) rf vesus <j>. Parameters are 
Ei = 1, E 2 = 0, E 3 = 1 , F = 1 , t = 0.4, Tl = 1.125, Tr = 1, fiL = —0.01, and fiR = 0. 



FIG. 10. Time-reversal symmetry breaking in a three-terminal triple QDs model for thermoelectrics, (a) TRB parameter r and 
(b) degree of coupling q as functions of the QDs energies Ei and E 2 . (c) The function 1 + q 2 r — q 2 — q 2 r 2 for the same values 
of energies as in panel (a), demonstrating that the quantum bound is satisfied. In (a)-(c) we used the following parameters: 
the QD connected to the probe terminal has energy E 3 = 2, cj> = 7r/2, F = 0.5, and t = —0.2 (energy unit is kpT = 1). (d) 
The example illustrates the appearance of r —>• 00 at a certain flux (0) value, when the parameters are chosen as E\ = 0.2, 
E 2 = —0.2, E 3 = 1, r = 0.5, and t = —1. 


and obtain the LDF of efficiency fluctuations as [2] 

g(n) = - minJ 2 t(-ii)C, - Vr)/T, A u A 2 ) 


(42) 


for any given A\ and A 2 . 

In Fig.^a) we plot the LDF J(rf) as obtained from the Gaussian approximation [Eq. (I24dl) ] and that from the full 
counting statistics method. The two functions perfectly match. In Fig. HKb) we plot the least probable efficiency 77 * 
calculated from the Gaussian approximation and the full counting statistics method. Again, the results from the two 
methods agree well with each other. These results confirms the validity of our analysis in the main text based on the 
Gaussian approximation, a consequence of the fluctuation theorem. 


IV. CHARACTERIZATION OF THERMOELECTRIC TRANSPORT IN A TRIPLE-QDS SYSTEM 


We calculate linear transport coefficients for a three-terminal thermoelectric model when the average thermal and 
electrical currents flowing out of the probe terminal P are zero. These conditions lead to M = M' LL — M' LP M'pp M' PL . 
Here M’ is the transport matrix for the total three-terminal system, i.e., /' = M'A' where /' = (I' L ,I p ) and 
A ' = {A'l,A'p) with = (4e> J 7 h) and A^ = (A( /e , Al /h ) are the currents and affinities for terminals 7 = L,P. I e 
is the charge current, Ih is the heat current, and e.g. , A Le = {n L - /ir)/T r , A Lh = 1 /T r - 1/T l . Linear transport 
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coefficients are calculated from the expression 

M' lY =J^ [V - |S 77 'M)| 2 3 4 5 6 7 8 ] ( l J ) /o(w)[l - /oM], (43) 

where (f>) (7,7' = L,P) is the scattering matrix between terminals 7 and 7', (j) = 27r$/$o- $ and $0 are the 

magnetic flux in our triple-quantum-dots (QDs) system and the flux quantum, respectively. The Fermi distribution 
/o(w) = [exp(y) + l ] -1 corresponds to an equilibrium state with the chemical potential set at fi = 0. Onsager 
reciprocal symmetry originates from the symmetry 5 77 '(w, <f>) = 

The transmission function is obtained from the relation S(uj,(/)) = —1 + iTG r (oj). Here 1 is a 3 x 3 identity 
matrix and G r (oj) = [(w + iT)l — H q d J ” 1 is the retarded Green’s function of the QDs. The hybridization energy 
T = 27 \ v k\ 2 S(oj — £k ) is assumed to be a constant (independent of energy) for all three QDs. 

We calculate the transport coefficients and then determine the TRB parameter r and the degree of coupling q. In 
Fig. [TO] we plot these parameters against E\ and £ 2 - We further confirm the “quantum bound” on linear transport 
coefficients discovered by Brandner et a 1. for three-terminal TRB conductors [g], which using our parametrization 
casts into the form as 1 + q 2 r — q 2 — q 2 r 2 > 0. In Fig. [TOT cl we plot it and show that it is always greater than zero in 
our parameter region. To illustrate the divergence of r, we plot it in Fig. HOf dl as a functions of the magnetic flux q i. 
We find that r diverges when M21 goes to zero with a finite M 12 . This happen when (j) ~ 0.47T or (j> ~ 1.67T. 
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